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1.  Introduction 


Research  is  being  conducted  at  the  Harry  Diamond  Laboratories  (HDL)  on 
millimeter-wave  radar  target  signatures.  A  coherent,  frequency-agile,  two 
coordinate,  amplitude  comparison  95-GHz  monopulse  radar  was  designed 
and  built  at  HDL  for  this  purpose.  To  simulate  the  data  that  will  be  obtained 
with  this  radar,  a  computer  program  was  written  that  calculates  the  signal 
return  from  a  complex  array  of  point  reflectors.  This  program  can  simulate 
frequency  agility,  both  AM  and  FM  system  noise,  and  signature  data  gener¬ 
ated  with  a  rotating  target  array.  Section  2  of  this  report  describes  the  data 
simulation  program. 

Two  additional  programs  were  written  to  analyze  the  data  generated  by  the 
simulation  program.  These  programs  form  the  basis  of  some  of  the  analysis 
programs  that  will  be  required  for  processing  data  obtained  with  the  radar. 
The  First  program  uses  a  Fast  Fourier  Transform  to  obtain  the  down-range 
location  of  the  simulated  target  array,  and  the  sum  and  difference  patterns 
generated  by  the  monopulse  radar  to  obtain  the  cross-range  location.  The 
second  program  uses  frequency- stepped  data  generated  as  a  target  array  is 
rotated  on  a  turntable  to  generate  a  two-dimensional  radar-cross-section- 
density  image  of  the  target  array.  This  is  done  with  Inverse  Synthetic  Aper¬ 
ture  Radar  (ISAR)  techniques.  These  programs  are  described  in  section  3  of 
this  report. 

In  section  4  of  this  report  we  discuss  the  results  of  using  the  programs  of 
section  3  with  the  data  generate  d  by  the  program  in  section  2.  The  discussion 
includes  high-range-resolution,  monopulse,  and  ISAR  results. 


2.  Monopulse  Radar  Signals 


2.1  Monopulse  Signals  From  Reflector  Arrays 


Consider  a  four-beam  two-coordinate  amplitude  comparison  monopulse 
radar  illuminating  an  array  of  M  point  reflectors.  The  radar  is  assumed  to  be 
located  at  the  origin  of  a  coordinate  system  and  the  f'  reflector  is  located  at 
position  P  (x  ,y.,z )  (see  fig.  1).  If  the  phase  of  the  transmitted  signal  is  0o  and 
that  of  the  received  signal  is  0,  then  the  shift  in  phase  of  a  signal  propagating 
from  the  radar  to  the  reflector  and  back  is  given  by 

A0  =  0  -  0O  (1) 


and  equals 

A8  =  +  yf  +  zi}'a  .  (2) 

where  c  is  the  speed  of  light  and/t  is  the  frequency. 


5 


Figure  1.  Monopulse 
calculation  beam 
geometry. 


The  amplitude  of  the  received  signal  is  dependent  on  the  distance,  I  Pj, ,  of  the 
reflector  from  the  origin,  the  reflectivity,  p.,  of  the  reflector,  and  on  iis  position 
in  the  radar  beam.  We  assume  that  each  beam  of  the  monopulse  radar  has 
circular  symmetry  about  the  individual  beam  axes  and  that  in  the  far  field  of 
the  antenna  the  normalized  field  intensity  at  an  angle,  <J>,  from  this  axis  is  the 
known  antenna  pattern,  £(<£>).  If  B  is  a  vector  with  the  same  magnitude  as  P 
pointed  along  the  axis  of  the  i'h  beam  of  a  monopulse  radar,  then  the  angle 
between  the  i*  beam  axis  and  the  /*  reflector  is  given  by 


cos  (<by)  = 


Pj  Bj 

W 


(3) 


The  field  intensity  for  a  point  P.  is  then  calculated  from  the  angle  dv  for  each 
beam. 


In  figure  1,  the  monopulse  axis  is  taken  to  be  along  the  y-axis  and  the 
individual  beams.  S.’s,  are  numbered  1  to  4.  The  squint  angle,  as,  is  the  angle 
between  each  Bt  and  the  x-y  or  the  y-z  plane. 

The  electric  field  intensity  at  the  receiver  for  the  monopulse  sum  is  therefore 


<=i;=i  Kj) 


f  f  p mf 

~kh  n1 


cos 


(¥w) 


(4) 


+  i  sin 


47C 


=  Az  +  iQn. 


If  the  initial  phase,  90,of  the  transmitted  beam  is  chosen  to  be  zero,  is  the 
in-phase  and  is  the  quadrature-phase  portion  of  the  received  sum  field 
intensity.  Similarly,  the  azimuth  and  elevation  differences, andX^ ,  are 
given  by 
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Equations  (4),  (5),  and  (6)  form  the  basis  for  the  monopulse  signal  calcula¬ 
tions. 

The  normalized  electric  field  pattern  of  each  beam,  £(F),  is  taken  to  be  that 
for  a  cosine  aperture  distribution  [1]: 

£(q.)  _  n  [sin  (v  +  n/2)  (  sin  (y  -  n/2)j  (7) 

4  [  y  +  k/2  \|/  -  7t/2 

where  vy  =  (  a/X)  sin  <t>,  a  is  the  aperture  diameter,  and  X  is  the  wavelength. 
The  wavelength  used  for  the  calculation  was  3.16  mm  and  the  aperture 
diameter  was  15.75  cm.  The  squint  angle,  as,  was  set  at  7.44  mrad,  slightly  less 
than  the  value  for  a  maximurrf  monopulse  error  slope  [1]. 

2.2  Data  Simulation  Program 

The  data  simulation  program  uses  the  equations  generated  in  the  previous 
section  to  calculate  the  monopulse  radar  signal  returned  from  an  array  of  point 
reflectors.  The  main  program  is  composed  of  two  sections.  The  first  calculates 
the  signal  return  for  a  frequency-agile  monopulse  radar  illuminating  a 
stationary  target  array.  The  second  extends  this  calculation  to  that  for  a 
rotating  target  array.  Both  sections  also  allow  for  the  addition  of  amplitude  or 
phase  noise  to  the  ideal  signal  calculations  to  simulate  radar  data  with  a  finite 
signal-to-noise  ratio. 
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The  stepped-frequency  monopulse  signal  calculations  are  made  using  the 
program  listed  in  appendix  A.  The  program  requires  the  user  to  input  the 
positions  of  each  point  of  the  target  array  (geometric  parameters)  and  the 
reflectivity  of  each  point  (target  parameters).  The  user  also  enters  the  number 
of  frequency  steps,  N,  and  the  frequency  step  size.  The  program  calculates  the 
summations  given  by  equations  (4),  (5),  and  (6)  for  N  frequencies  beginning 
at  94.44  GHz.  The  values  calculated  for  the  in-phase  and  quadrature 
components,  /and  Q,  for  each  frequency  are  placed  in  a  complex  array  which 
is  then  stored  in  a  disc  file.  The  data  in  this  file  are  used  in  separate  programs 
for  analyses  that  will  be  discussed  in  the  following  sections  of  this  report. 

The  second  section  of  the  simulation  program  generates  monopulse  radar 
signals  for  an  array  of  point  isotropic  reflectors  on  a  turntable.  In  this  case,  the 
target  array  and  the  radar  line-of-sight  are  assumed  to  lie  in  the  plane  of  the 
turntable.  The  radar  line-of-sight  also  is  taken  to  pass  through  the  axis  of 
rotation,  which  is  perpendicular  to  the  plane  of  the  turntable.  With  the  data 
generated  by  this  program,  ISAR  techniques  can  be  used  to  calculate  high- 
resolution  images  of  the  target  in  both  down-range  and  cross-range  configu¬ 
rations. 

The  program  starts  the  calculation  by  recording  the  amplitude  and  phase 
history  of  a  point  reflector  as  it  moves  on  a  circle  about  the  center  of  rotation. 
The  value  of  I  P\  in  equation  (4)  is  computed  in  the  following  manner.  The 
range  from  the  radar  to  the  center  of  rotation  of  the  target  array  is  defined  by 
the  vector  yg,  the  radial  distance  of  the  j,h  point  to  the  center  of  rotation  is  taken 
to  be  r ,  and  the  initial  angle  of  each  point  reflector  is  w..  The  angle  w0.  is  the 
angle  between  r  and  a  line  perpendicular  toyc  and  passing  through  the  center 
of  rotation.  The  simulation  program  computes  I  P\  for  each  point  reflector  and 
from  it  the  radar  signal  (eq  (4),  (5),  and  (6))  for  N  frequency  steps.  Each  w 
is  then  incremented  by  an  angle.  Aw,  where  Aw  is  defined  to  be  negative  for 
clockwise  rotation,  and  the  calculation  is  repeated.  The  resulting  P.  following 
m  angular  increments  is 

| Pj\  =  [rf  +yo2  +  2  rjy0  sin  (wOJ  +  mAw)]1/2  .  (8) 

The  calculation  is  then  repeated  for  as  many  aspect  angles  as  is  desired  to 
generate  the  ISAR  image.  The  simulation  consisted  of  64  sets  of  different 
aspect  angles  separated  by  Aw  with  64  frequency  steps  of  10  MHz  for  each 
aspect  angle.  The  calculated  /,  Q  values  for  each  frequency  and  each  position 
were  stored  into  a  two-dimensional  complex  array.  This  array  was  then  stored 
in  a  file  to  be  processed  by  another  program,  which  will  be  discussed  in  section 
3.3. 

Once  the  basic  signal  calculations  are  made,  another  section  of  the  program 
allows  the  addition  of  different  types  of  noise  to  the  calculated  monopulse 


sum  and  difference  signals  to  simulate  actual  measurement  conditions.  For 
a  frequency  fk  =  /,  the  equations  for  the  sum  and  difference  /,  Q  values  (eq  (4), 
(5),  and  (6))  were  modified  and  expressed  in  the  form  shown  below  to  add 
various  types  of  noise  components: 


Iz  =  (A+  A')  cos  (0  +  O')  +  D'  cos  03 

(9) 

Qz  =  (A  +  A')  sin  (0  +  0')  +  D'  sin  03 

(10) 

f\4Z  =  6  cos  0i  +  B'  cos  04 

(ID 

Qmz  =  B  sin  0i  +  B'  sin  ©4 

(12) 

/a el  =  C  COS  02  +  C'cos  05 

(13) 

Qsel  =  C  sin  02  +  C'sin  ©5 

(14) 

where  A,  B,  C,  0,  0,,  and  02  are  the  sum,  azimuth  difference,  and  elevation 
difference  signal  amplitudes  and  phases,  respectively.  Pulse-to-pulse  ampli¬ 
tude  and  phase  noise  are  represented  by  the  Gaussian-distributed  variables  A ' 
and  0',  which  have  been  added  to  only  the  sum  channel  signals  in  our 
simulation.  For  the  difference  channel  signals,  the  dominant  noise  is  system 
noise,  which  is  set  by  the  system  noise  figure.  This  is  so  because  the  difference 
channel  signals  are  normally  small  for  on-axis  targets,  producing  a  much 
smaller  signal-to-noise  ratio  in  the  difference  channels.  In  addition,  both 
pulse-to-pulse  amplitude  and  phase  noise  will  affect  all  the  channels  equally, 
so  that  the  normalized  difference  signal,  which  determines  off-axis  distance, 
will  not  be  affected  by  the  addition  of  pulse-to-pulse  noise.  Noise  contribu¬ 
tions  due  to  poor  system  signal-to-noise  ratio  and  random  phase  noise  are 
represented  by  the  random  noise  amplitude  components  D',B\  and  C '  and  the 
random  phase  components  03,  04,  and  0y 

The  different  types  of  noise  sources  were  compared  by  running  simulations 
with  each  noise  separately  and  in  one  case  with  all  types  of  noise  included. 
Simulations  were  run  with  AM  noise  for  pulse-to-pulse  signal-to-noise  ratios 
of  30,  20,  and  10  dB.  System  signal-to-noise  ratios  of  40,  30,  20,  and  10  dB 
were  used  to  analyze  the  effect  of  poor  signal-to-noise  ratio  on  system  per¬ 
formance.  In  addition,  pulse-to-pulse  phase  errors  of  5,  10,  and  20  deg  were 
added  to  the  signal  in  another  set  of  simulations.  The  final  simulation  was  run 
using  a  system  signal-to-noise  ratio  of  20  dB,  a  pulse-to-pulse  signal-to-  AM- 
noise  ratio  of  20  dB  and  a  pulse-to-pulse  phase  error  of  5°. 

The  data  generated  in  these  simulations  are  used  in  data  processing  programs 
that  generate  high-range-resolution  images  and  two-dimensional  radar  IS  AR 
images  of  the  target  arrays.  In  addition  the  data  are  used  in  calculating  the  off- 
axis  monopulse  response  for  the  simulated  radar.  The  results  of  these 
simulations  will  be  discussed  in  section  4  of  this  report. 


3.  Data  Processing  Techniques 


3.1  Fourier  Transform  Processing 

The  basic  techniques  for  processing  the  frequency-stepped  monopulse  radar 
data  rely  on  the  use  of  the  Fast  Fourier  Transform  (FFT)  algorithm  to  generate 
high  resolution  range  gates.  For  evaluating  turntable  data,  a  two-dimensional 
FFT  can  be  used  to  generate  high  resolution  two-dimensional  ISAR  images 
of  the  target.  The  applications  of  both  of  these  processing  techniques  to  the 
simulated  radar  data  will  be  discussed  in  this  section. 


The  FFT  is  a  useful  tool  for  performing  both  time-domain  and  frequency 
domain  signal  analysis.  The  most  common  use  of  the  FFT  is  to  transform  time 
domain  signals  into  their  corresponding  spectral  components.  This  type  of 
transform  generates  a  function  with  amplitude  and  phase  values  over  a 
defined  frequency  range.  The  basic  Fourier  transform  pair  equations  are  given 
by  the  following  integrals  [2]: 


■I 


=  I  X(f)  exp  (pLnfi)  df 


« — 


X{f)  =  jit)  exp  (~J2vft)  dt  , 


05) 


(16) 


where  the  first  integral  defines  the  inverse  Fourier  transform  and  the  second 
integral  the  direct  Fourier  transform.  For  sampled  data,  the  discrete  Fourier 
transform  (DFT)  is  used.  These  transforms  are  defined  by  the  following 
summations  [21: 


N- 1 

x(nAt)  =  A/^r  X^/:A/)exp (jlnkAfnAt)  ,  and  (17) 

k=0 

N- 1 

X J,kAf)  =  At  ^  x(nAt)e\p  (-jlnkAfnAt)  ,  (18) 

n= o 


where 


N 

At 

A/ 

n 

k 

.r(flAf) 

X(kAf) 


number  of  samples, 

time  between  samples, 

sample  interval  in  the  frequency  domain, 

the  time  sample  index  n  =  0,  1 ,  2...  N- 1 , 

the  index  for  discrete  frequency  components,/:  =  0, 1,2...  AM . 

discrete  set  of  time  samples  that  defines  the  waveform, 

the  set  of  Fourier  coefficients  obtained  by  the  DFT  of  x{nf). 
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These  equations  form  the  basis  for  computer  calculations  of  the  Fourier 
transform.  Many  computers  have  a  built-in  FFT  algorithm  that  is  much  faster 
than  a  direct  calculation  of  the  summations.  This  analysis  is  not  intended  to 
be  a  dissertation  on  the  Fourier  transform,  but  we  will  use  the  properties  of  the 
Fourier  transform.  These  properties  are  discussed  in  greater  detail 
elsewhere  [2], 

Since  the  Fourier  transform  has  an  inverse  transform,  it  is  apparent  that  one 
can  transform  from  the  frequency  domain  to  the  time  domain  as  well  as  from 
the  time  domain  to  the  frequency  domain.  With  a  frequency-agile  coherent 
system,  the  data  represent  information  in  the  frequency  domain,  and  when  the 
FFT  is  used,  the  data  are  transformed  into  the  time  domain.  Using  this 
technique,  high-resolution  range  gates  are  generated  from  wide-bandwidth 
stcpped-frequency  signals. 

In  transforming  the  data  from  the  frequency  domain  to  the  time  domain,  a 
rectangular  function  in  the  frequency  domain  is  transformed  into  a  sin  x/x 
distribution  in  the  time  domain  in  a  similar  way  as  a  rectangular  pulse  in  the 
time  domain  transforms  into  a  sin  x/x  distribution  in  the  frequency  domain.  To 
demonstrate  this  effect,  a  computer  simulation  was  run  to  calculate  the  signal 
return  for  two  targets  located  at  93.685  and  101.297  m,  each  having  a  radar 
cross  section  of  1  m2.  The  simulation  calculated  the  signal  return  for  64 
frequency  steps  of  10  MHz  beginning  at  94.44  GHz.  The  inverse  DFT  was 
used  to  generate  the  data  shown  in  figure  2(a).  The  plot  gives  the  calculated 
amplitude  of  each  scatterer  as  a  function  of  the  point  number  of  the  DFT. 

The  resolution  cell  width  is  given  by  [3] 

AR  =  cflNkf  ,  (19) 

where  c  is  the  speed  of  light,  N  is  the  number  of  frequency  steps,  and  Afis  the 
frequency  step  size.  For  64  steps  of  10  MHz  the  resolution  cell  width  or  range 
resolution  is  0.234  m.  The  total  range  coverage  of  the  DFT  is  given  by  [3] 

R  =  cfltxf  .  (20) 

For  the  1 0-MHz  step  size,  the  total  range  coverage  is  1 4.99  m.  This  means  that 
the  target  range  is  broken  up  into  range  intervals  approximately  1 5  m  in  extent 
that  overlap  every  1 5  m;  that  is,  a  target  at  1 0  m  will  image  at  the  same  position 
as  targets  at  15 n+  10  m,  where  n  is  an  integer.  Therefore,  this  technique  does 
not  give  absolute  range,  but  only  relative  range  within  the  14.99  ni  of 
coverage. 

It  is  apparent  from  fig.  2(a)  that  the  two  points  are  not  being  imaged  in  the 
same  manner.  This  occurred  because  the  first  or  closest  point  was  placed  in 
the  center  of  a  range  bin,  whereas  the  second  point  was  placed  between  two 
range  bins.  For  the  point  in  the  center  of  the  bin,  the  peak  amplitude  occurs 
at  the  center  of  a  range  cell,  and  the  samples  near  the  peak  occur  at  the  nulls 
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Figure  2.  Range 
imaging  using  64- 
point  FFT  with  a 
rectangular  window 
md  imaging  with 
padded  128-point 
FFT. 
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of  the  sin  x/x  distribution.  For  the  other  point  the  amplitude  is  reduced  by  about 
4  dB,  and  the  samples  around  the  peak  now  occur  at  the  maxima  of  the  side 
lobes  of  the  sin  x/x  distribution. 

This  amplitude  variation  is  caused  by  the  location  of  the  target  with  respect 
to  the  range  gates  generated  by  the  FFT.  For  a  target  located  at  the  center  of 
a  gate,  the  amplitude  response  is  maximum.  For  a  target  located  in  between 
gates,  the  amplitude  response  is  determined  by  both  gates.  To  reduce  the 
amplitude  variation,  overlapping  range  gates  with  gates  lying  in  between  the 
original  gates  need  to  be  generated.  This  is  simply  accomplished  by  extending 
the  original  record  length  with  a  set  of  samples  that  are  identically  zero.  For 
our  purposes,  32  zeros  are  added  at  the  beginning  and  the  end  of  the  data  set. 
This  process  is  referred  to  as  padding,  and  in  this  case  a  128-point  DFT  is 
performed.  The  results  of  this  procedure  are  shown  in  figure  2(b).  The 
amplitude  of  both  points  is  now  the  same,  but  for  a  point  located  between 
them,  there  would  still  be  an  amplitude  reduction.  There  is  also  an  overall 
reduction  in  amplitude  of  6  dB.  The  peaks  and  nulls  of  the  sinx/x  distribution 
are  now  easily  distinguishable.  To  further  enhance  the  response  and  reduce 
the  amplitude  variation,  more  zeros  could  be  added,  but  the  computation  time 
would  be  greatly  increased.  For  all  subsequent  data  analyses  in  this  report,  64 
zeros  were  added  to  the  data  before  processing. 

The  images  shown  in  figure  2(b)  can  be  improved  further  by  changing  the 
window  function  used  in  processing  the  data.  For  the  initial  data  processing, 
a  rectangular  window  function  was  used,  which  generated  sidelobes  begin- 
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ning  at  13  dB  down  from  the  peak.  In  generating  range  images  of  complex 
targets,  this  type  of  sidelobe  performance  would  greatly  restrict  the  usefulness 
of  the  images.  There  are  a  variety  of  other  potential  window  functions;  some 
of  the  most  common  are  listed  in  table  1.  The  table  includes  the  sidelobe 
performance,  corresponding  reduction  in  peak  level,  and  3-dB  mainlobe 
width  of  the  window  functions  where  B  =  1  /T.  For  this  application,  the 
Hamming  window  provides  the  best  sidelobe  performance.  Use  of  this 
window  function  is  demonstrated  in  figure  3  for  the  same  simulation  points 
as  in  figure  2.  The  FFT  point  number  for  this  Figure  and  all  subsequent  Figures 
has  been  replaced  by  relative  range,  since  the  point  number  is  proportional  to 
the  range  within  the  1 5-m  window.  Easily  discernible  differences  are  a  reduc¬ 
tion  in  the  First  sidelobe  level  of  29  dB,  a  reduction  in  peak  level  of  5.35  dB, 
and  an  increase  in  the  3-dB  width  of  the  main  response. 

In  the  remainder  of  this  report,  the  Hamming  window  will  be  used  in 
conjunction  with  padding  to  improve  the  down-range  imaging  properties  of 
the  data  analysis.  The  total  reduction  in  peak  response  due  to  the  windowing 
and  padding  is  11.35  dB,  and  the  down-range  images  will  reflect  this 
reduction  in  amplitude. 

3.2  Monopulse  Signal  Processing 

The  ideal  monopulse  radar  calculates  the  target  off-axis  angle  error  by  a 
comparison  between  the  difference  channel  signal  and  the  sum  channel 
amplitude  at  each  range  gate.  When  a  complex  target  with  range  extent  less 
than  the  range  cell  determined  by  the  transmitted  pulsewidth  is  illuminated 
with  a  frequency-agile  monopulse  radar,  the  target  can  be  resolved  in  range 
with  higher  resolution  by  using  the  DFT  algorithm.  The  absolute  range  is  then 
determined  by  the  time  delay  to  the  sampling  window,  and  the  relative  range 


Table  I.  Selected  Data 
Window  Types 

Showing  Performance 

Window  type 

Maximum 
sidelobe  (dB) 

Peak  ampli¬ 
tude 

Mainlobe 
3-dB  width 

Specification 

Rectangular 

-13.2 

T 

0.86  B 

Triangular 

-26.7 

0.5T 

1.27  B 

Cosine  (Hanning) 

-31.6 

0.5T 

1.39  B 

Hamming 

-41.9 

0.54T 

1.26  B 

Figure  3.  Imaging 
with  padded  128- 
point  FFr  using 
Hamming  window 
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within  that  window  is  further  resolved  using  the  DFT.  In  addition,  tne 
difference-channel  data  can  be  processed  using  the  DFT  to  provide  a  down- 
range  profile  of  both  azimuth  and  elevation  error  signals.  By  applying  a  cross¬ 
range  calibration  to  these  error  signals,  we  can  determine  the  cross-range 
position  of  the  scattering  center. 

To  generate  a  monopulse  calibration  curve  for  the  off-axis  angle  error,  a 
monopulse  simulation  was  run  using  a  single-point  scatterer  as  the  target.  The 
target  was  positioned  at  a  range  of  100  m,  at  the  center  of  the  elevation  axis, 
and  translated  through  an  azimuth  angle  of  -20  to  +20  mrad  from  the 
monopulse  axis.  The  azimuth  error  signal  was  calculated  along  with  the  sum 
channel  signal.  An  off-axis  calibration  curve  was  generated  by  normalization 
of  the  difference  channel  signal  to  the  sum  channel  signal.  The  results  of  this 
calculation  are  shown  in  figure  4.  Since  the  monopulse  beam  pattern  was 
assumed  to  be  symmetric  in  both  elevation  and  azimuth,  the  normalized 
difference  channel  calibration  curve  can  be  used  to  calculate  the  angle  error 
for  both  axes.  The  results  of  applying  this  technique  to  the  analysis  of  a 
complex  target  are  discussed  in  section  4.2  of  this  report. 


3.3  ISAR  Processing 


The  basic  imaging  technique  used  in  generating  a  two-dimensional  target 
image  from  coherent  frequency-agile  rotating  target  data  is  very  similar  to 
that  employed  in  generating  a  synthetic  aperture  radar  (SAR)  image.  In  the 
latter  method,  the  radar  is  moved  around  or  past  the  target,  and  the  return 
signal  amplitude  and  phase  are  recorded  as  a  function  of  radar  position.  The 
movement  of  the  radar  synthesizes  a  larger  antenna  aperture,  which  improves 
the  radar  resolution.  In  the  ISAR  method,  the  target  is  rotated  and  the  radar 
is  held  fixed  while  the  amplitude  and  phase  history  of  the  target  are  recorded. 
To  generate  the  image,  a  two-dimensional  FFT  is  used.  The  FFT  is  performed 
on  the  frequency-stepped  data  to  provide  a  down-range  amplitude  and  phase 
history  of  the  individual  scatterers  as  they  are  rotated  through  a  fixed  angle. 
Included  in  the  phase  history  is  the  sampled  Doppler  frequency  of  each 
scatterer  due  to  its  relative  motion  with  respect  to  the  radar.  The  second  FFT 
is  used  to  analyze  the  Doppler  frequencies,  both  positive  and  negative,  that 


Figure  4.  Monopulse 
sum  and  normalized 
azimuth  difference 
pattern. 
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exist  in  the  phase  history.  These  Doppler  frequencies  can  be  related  directly 
to  the  off-axis  distance  of  the  individual  scatterers.  The  details  and  applica¬ 
tions  of  high  resolution  imaging  techniques  have  been  reported  elsewhere 
[3-6].  We  will  use  these  results  and  discuss  the  application  of  these  methods 
in  analyzing  the  simulated  and  actual  radar  data. 

To  generate  the  IS  AR  image,  a  set  of  frequency-stepped  data  are  generated  for 
many  different  aspect  angles.  For  our  simulations,  64  different  frequencies 
make  up  the  frequency  stepped  data,  and  64  different  aspect  angles  are  used 
for  generating  the  cross-range  resolution.  The  down-range  resolution  is  deter¬ 
mined  by  the  frequency-stepped  data  as  described  earlier,  whereas  the  cross¬ 
range  resolution  is  determined  by  the  total  angle  of  observation.  The  cross¬ 
range  resolution  is  given  by  [3] 

Ax  =  c/lFMAw  ,  (21) 

where  c  is  the  speed  of  light,  F  is  the  average  frequency  during  a  frequency 
ramp,  M  is  the  number  of  different  aspect  angles  observed,  and  Aw  is  the 
angular  increment  between  successive  frequency  ramps.  The  total  cross¬ 
range  distance  imaged  is  given  by  [3] 

x  =  c/2FAw  .  (22) 

There  are  some  restrictions  that  must  be  adhered  to  during  the  generation  of 
the  data  if  the  image  quality  is  to  be  preserved.  The  rotation  rate  has  an  effect 
on  the  quality  of  both  the  down-range  and  cross-range  images.  If  the  rotation 
rate  is  too  high,  there  will  be  significant  down-range  and  cross-range  image 
degradation,  as  discussed  next. 

The  down-range  image  is  generated  from  the  phase  history  of  the  target  as  it 
moves  during  the  frequency  stepping.  For  a  stationary  target,  the  target  is 
imaged  accurately,  but  if  there  is  target  motion,  the  Doppler  frequency  can 
distort  the  down-range  position.  For  accurate  imaging,  the  Doppler  shift  from 
target  motion  must  be  kept  below  the  effective  sampling  frequency  of  the 
frequency-stepped  waveform.  The  Doppler  shift  from  target  motion  is  given 
by 

Fd  =  2ra/k  ,  (23) 

where  r  is  the  target  radius,  0)  is  the  rotation  rate  in  radians  per  second,  and 
X  is  the  wavelength.  The  effective  sampling  frequency  of  the  frequency- 
stepped  waveform  is  given  by 

F,  =  PRF/N  ,  (24) 

where  PRF  is  the  pulse  repetition  rate  and  N  is  the  number  of  steps  required 
to  generate  a  down-range  image.  To  maintain  a  reasonable  image  quality,  F d 
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must  be  kept  less  than  one-eighth  of  Fr  This  means  that  the  rotation  rate  is 
limited  by 


U<PRF)in6rN  .  (25) 

For  a  radar  operating  at  a  wavelength  of  3.16  mm,  using  64  frequency  steps, 
and  a  repetition  rate  of  10  kHz,  and  with  a  maximum  target  radius  of  5  m,  the 
maximum  rotation  rate  is  equal  to  0.0062  rad/s  or  one  complete  revolution  of 
the  target  every  17  min.  In  the  rotating  target  simulation  reported  here,  the 
frequency-stepped  data  were  generated  while  the  target  was  stationary,  so 
there  is  no  apparent  Doppler  frequency  included  in  the  stepped  frequency 
data.  The  target  aspect  angle  was  incremented  by  0.0088 1 0  between  ramps  to 
generate  the  cross-range  resolution.  The  Doppler  is  due  to  the  motion  between 
frequency  ramps,  and  it  is  contained  in  the  aspect  angle  data.  The  cross-range 
resolution  is  0.16  m. 

The  ISAR  data  are  processed  in  the  following  manner.  The  amplitude  and 
phase  are  calculated  for  each  frequency  as  the  frequency  is  stepped,  and  the 
64  data  points  are  stored  in  a  complex  array.  This  procedure  is  repeated  for 
each  aspect  angle,  thereby  generating  a  64  x  64  complex  array.  This  complex 
array  is  then  windowed  in  both  dimensions  with  a  Hamming  window  being 
used  to  reduce  sidelobes.  The  data  are  then  padded  with  32  zeros  at  the 
beginning  and  end  of  each  frequency  ramp,  and  further  padded  with  32  ramps 
of  zeros  at  the  beginning  and  end  of  the  aspect  angle  data.  This  generates  a 
complex  array  with  1 28  x  128  elements.  A  DFT  algorithm  is  applied  to  the  128 
sets  of  frequency-stepped  data,  and  the  resulting  array  is  then  transformed  by 
aspect  angle.  The  array  must  be  inverted  to  order  the  Doppler  frequency  data, 
which  correspond  to  the  cross-range  coordinate,  to  produce  the  proper  image. 
The  resulting  amplitude  data  are  corrected  by  adding  22.8  dB  to  compensate 
for  the  reduction  in  signal  amplitude  due  to  windowing  and  padding.  The 
array  is  plotted  in  down-range  and  cross-range  distance  to  give  a  two- 
dimensional  radar-cross-section  image  of  the  target.  The  results  of  this 
simulation  will  be  discussed  in  the  next  section. 


4.  Data  Simulation  Results 

4.1  High-Range-Resolution  Results 

The  range  imaging  capability  of  FFT  processing  was  shown  in  Figure  3.  This 
section  will  discuss  the  effects  of  different  types  of  noise  sources  on  the 
quality  of  images  that  can  be  generated  with  FFT  processing.  A  baseline 
simulation  was  run  with  a  point  reflector  (with  radar  cross  section  of  1  m2) 
placed  97.433  m  from  the  radar.  The  data  were  generated  as  discussed  in 
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section  3. 1 .  Figure  5  shows  the  results  of  the  simulation.  The  amplitude  of  the 
peak  is  reduced  by  1 1 .35  dB  because  of  windowing  and  padding,  with  a  range 
response  indicative  of  the  Hamming  window  function.  The  calculations  were 
performed  with  double  precision  arithmetic  to  minimize  processing  errors 
due  to  truncation,  and  there  were  no  noise  components  added  to  the  signal. 
This  figure  will  be  the  baseline  forcomparing  the  effects  of  the  different  types 
of  noise  to  the  ideal  case. 

The  first  type  of  noise  considered  is  that  from  pulse-to-pulse  amplitude 
fluctuations,  which  can  be  considered  to  be  pure  AM  noise.  The  results  of 
adding  this  kind  of  noise  to  the  calculated  data  are  shown  in  figures  6,  7,  and 
8.  The  AM  signal-to-noise  ratios  are  30, 20,  and  10  dB,  respectively;  this  ratio 
was  vafied  by  changing  the  noise  contribution  while  keeping  the  signal  level 
constant.  The  figures  indicate  an  increase  in  the  baseline  noise  level  as  the 
signal-to-noise  ratio  is  changed  from  30  to  20  to  10  dB.  This  increase  is 
proportional  to  the  change  in  signal-to-noise  ratio,  but  there  is  no  apparent 
change  in  the  peak  amplitude  response.  For  the  10-dB  fluctuation  case,  the 
amplitude  of  the  peak  response  still  exceeds  the  average  noise  level  by  about 
28  dB.  This  demonstrates  the  signal-to-noise  ratio  improvement  that  can  be 
generated  by  coherently  processing  64  different  samples  of  the  same  signal 
plus  Gaussian  distributed  noise.  Figure  6  shows  that  a  fluctuation  ratio  of  30 
dB  does  not  produce  significant  degradation  from  the  noiseless  case,  since  the 
noise  level  is  almost  equivalent  to  the  maximum  sidelobe  level  generated  with 
the  Hamming  window. 


Figure  5.  Baseline 
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Figure  6.  Effect  of  -30 
dB  pulse-to-pulse 
amplitude  noise  below 
signal  on  128-point 
FFT  single-point 
reflector  range  image. 
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Figure  7.  Effect  of 
-20  dB  pulse-to-pulse 
amplitude  noise 
below  signal  on  128- 
point  FFT  single¬ 
point  reflector  range 
image. 
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Figure  8.  Effect  of  -10 
dB  pulse-to-pulse 
amplitude  noise  below 
signal  on  128-point 
FFT  single-point 
reflector  range  image. 


Relative  Range  (m) 

Another  interesting  feature  of  these  figures  is  the  symmetry  of  the  baseline 
noise  around  the  peak  response.  For  pure  AM  noise,  the  detected  phase  is  not 
affected  by  the  noise  modulation.  Since  the  detected  phase  determines  the 
location  of  the  target  return,  the  noise  will  simulate  a  fluctuating  target  at  the 
real  target  location.  The  resulting  distribution  will  be  symmetrical  in  time 
about  the  target  location.  Pure  AM  noise  will  be  only  one  component  of  the 
total  noise  in  a  real  radar  system,  but  this  simulation  is  useful  in  determining 
how  much  amplitude  noise  is  acceptable. 

The  second  type  of  noise  to  be  modeled  is  system  noise.  As  described  earlier, 
system  noise  was  simulated  by  adding  a  random  amplitude  fluctuation  with 
a  random  phase  component  to  both  the/,  (/-channels.  System  noise  will  affect 
all  data  channels  randomly  and  will  be  the  dominant  noise  factor  for  small 
signal  levels.  Sample  simulations  were  run  for  signal-to-noise  ratios  of 40, 30, 
20,  and  10  dB.  The  results  are  shown  in  figures  9  through  12,  respectively. 
These  figures  closely  match  those  done  for  pulse-to-pulse  amplitude  noise, 
except  that  the  noise  distribution  is  no  longer  symmetrical  about  the  point 
location.  This  effect  is  caused  by  the  variation  in  the  total  signal  phase  created 
by  adding  both  a  random  amplitude  and  random  phase  noise  modulation.  A 
with  the  pulse-to-pulse  case,  the  noise  does  not  cause  much  image  degrada¬ 
tion  even  at  lOdB,  but  there  is  a  significant  loss  of  dynamic  range  in  the  image 
when  the  noise  level  exceeds  20  dB. 
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The  final  type  of  noise  which  was  modeled  separately  was  that  caused  by 
random  phase  errors.  These  can  be  caused  by  loss  of  phase  lock,  by  target  or 
radar  motion  during  data  acquisition,  or  by  phase  errors  in  the  /,  Q  detection 


Figure  9.  Effect  of  40- 
dB  system  signal-to- 
noise  ratio  on  128- 
point  FFT  single-point 
reflector  range  image. 


Figure  10.  Effect  of  30- 
dB  system  signal-to- 
noise  ratio  on  128-point 
FFT  single-point 
reflector  range  image. 


Figure  1 1.  Effect  of  20- 
dB  system  signal-to- 
noise  ratio  on  128-point 
FFT  single-point 
reflector  range  image. 
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Figure  12.  Effect  of  10- 
dB  system  signal-to- 
noise  ratio  on  128-point 
FFT  single-point 
reflector  range  image. 


system.  The  first  two  types  of  phase  error  were  modeled  by  adding  a  random 
phase  to  the  phase  calculated  using  the  Pt  defined  by  equation  (8).  This  phase 
error  is  common  to  all  the  data  channels,  so  both  the  sum  and  difference 
channels  are  affected  equally.  Simulations  were  run  with  random  phase  errors 
of  5, 10,  and  20°;  the  results  are  shown  in  figures  13, 14,  and  15.  The  resulting 
down-range  images  are  very  similar  to  those  generated  by  system  noise.  In 
fact,  the  phase-error  plot  for  a  5°  phase  error  is  almost  identical  to  the  system 
noise  plot  for  a  20-dB  signal-to-noise  ratio. 
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Figure  13.  Effect  of  5° 
pulse-to-pulse  phase 
variation  on  128-point 
FFT  single-point 
reflector  range  image. 


The  phase  noise  caused  by  errors  in  the  I,  Q  detection  system  was  not  modeled 
separately  because  its  effect  is  identical  to  that  caused  by  the  system  noise. 
However,  a  simulation  was  run  using  a  combination  of  the  AM  noise,  system 
noise,  and  random  phase  noise.  The  components  used  for  this  simulation  were 
20  dB  for  the  AM  noise  contribution,  20  dB  for  the  system  noise  contribution, 
and  5°  for  the  phase  noise  contribution.  The  results  of  this  simulation  are 
shown  in  figure  16.  The  combination  of  all  the  noise  sources  results  in  a 
somewhat  higher  background  noise  level  than  any  of  the  individual  cases. 

Table  2  shows  the  results  of  averaging  the  background  noise  levels  for  all  the 
noise  simulations.  The  peak  signal  level  for  each  simulation  was  -11.35 
dBsm,  so  that  the  resulting  signal-to-noise  level  can  be  calculated  by  dividing 
by  the  average  noise  levels.  The  results  of  this  calculation  are  also  shown  in 
table  2.  The  signal-to-noise  ratio  improvement  factor  is  about  1 8  dB  for  all  the 
cases  shown  in  table  2.  For  the  combination  case,  the  three  different  types  of 
noise  add  up  to  produce  three  times  as  much  noise  as  the  individual  compo¬ 
nents.  This  is  demonstrated  by  an  increase  of  the  average  noise  for  this  case 
of  approximately  5  dB  over  the  individual  cases. 


Figure  14.  Effect  of  10° 
puLse-to-pulse  phase 
variation  on  128-point 
FFT  single-point 
reflector  range  image. 
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Figure  15.  Effect  of  20° 
pulse-to-pulse  phase 
variation  on  128-point 
F  FT  single-point 
reflector  range  image. 
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Figure  16.  Effect  of  a 
combination  of  20-dB 
system  signal-to-noise 
ratio,  -20-dB  pulse- 
to-pulse  amplitude 
noise,  and  5°  pulse-to- 
pulse  phase  variation 
on  128-point  FFT 
single-point  range 
image. 
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Table  2.  Effect  of 
Different  Types  of 
Noise  on  High- 

Noise  type 

Noise  amplitude 
(relative  to  signal) 

Output  noise 
(Avg.  dBsm) 

Output 

S/N  (dB) 

Resolution-Range 

System 

-40  dB 

-70.6 

59 

FFT  Imaging 

-30  dB 

-58.4 

47 

-20  dB 

-49.3 

38 

-10  dB 

-40.8 

29 

Pulse-to-pulse 

-30  dB 

-60.5 

49 

-20  dB 

-50.7 

39 

-lOdB 

-39.8 

28 

Phase 

5° 

-50.2 

39 

10° 

-44.4 

33 

20° 

-38.8 

27 

Combination 

-20dB,-20  dB, 

-44.1 

33 

_ 5^ _ 

4.2  Monopulse  Results 

An  array  of  five  isotropic  point  reflectors  was  used  for  the  monopulse 
simulations.  The  reflectors,  each  with  a  radar  cross  section  of  1  m2,  were 
located  at  the  following  down-range  and  cross-range  positions:  point  No.  1  at 
92.93  m  and  0.0  m;  point  No.  2  at  94.93  m  and  + 1 .5  m;  point  No.  3  at  97.43 
m  and  +2.0  m;  point  No.  4  at  97.43  m  and  -2.0  m;  and  point  No.  5  at  101 .93 
m  and  -1 .0  m.  The  positive  cross-range  positions  indicate  the  left  side  of  the 
azimuth  beam  axis,  and  the  negative  positions  are  for  points  to  the  right  side 
of  the  azimuth  axis.  The  elevations  of  the  reflectors  were  all  the  same, 
although  the  simulation  could  have  included  variations  of  the  elevation 
positions  as  well.  This  target  array  also  was  used  in  the  IS  AR  calculations  that 
will  be  discussed  in  the  next  section. 

The  monopulse  off-axis  error  can  be  calculated  for  each  transmit  pulse  or  by 
integrating  many  pulses.  For  a  long-range  monopulse  tracking  radar,  either  of 
these  methods  can  be  used  to  calculate  the  off-axis  error.  For  the  application 
discussed  here,  the  high-resolution  imaging  capability  of  a  frequency-stepped 
radar  can  be  used  to  further  resolve  the  target  in  cross-range  as  well  as  in 
down-range,  but  only  if  there  is  only  one  target  in  each  range  bin. 
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Figure  17.  Effect  of  40- 
dB  system  signal-to- 
noise  ratio  on 
determining  off-axis 
target  position  using 
monopulse  calibration 
on  high-range- 
resolution  target  image 
data. 


Monopulse  returns  were  calculated  with  system  signal-to-noise  ratios  of  40, 
30,  and  20  dB.  The  results  of  processing  these  simulations  using  the  FFT 
algorithms  described  earlier  are  shown  in  figures  17, 18,  and  19.  The  top  parts 
of  the  figures  show  the  down-range  image  generated  by  FFT  processing  the 
sum  pattern  response  of  the  radar.  This  signal  was  processed  in  the  same 
manner  as  that  of  figure  5.  The  lower  parts  of  the  figures  show  the  cross-range 
response  generated  by  processing  the  difference  channel  signal  in  the  follow¬ 
ing  manner.  First,  the  difference  channel  signal  was  processed  using  the  FFT 
to  provide  a  high-resolution  down-range  difference  channel  response.  The 
sum  channel  signal  was  thresholded  at  —40  dBsm  to  determine  the  location  of 
the  different  scatterers.  The  difference  channel  signal  was  then  compared  to 
the  sum  channel  signal  point  by  point  to  determine  the  normalized  off-axis 
signal  level.  The  calibration  shown  in  figure  4  was  used  to  compute  the  off- 
axis  position  of  the  target  in  milliradians. 

The  calculated  positions  of  the  targets  are  in  excellent  agreement  with  the 
input  positions.  The  target  that  was  positioned  on-axis  does  not  show  up  in  the 
off-axis  response,  and  the  pair  of  points  positioned  symmetrically  about  the 
axis  also  appear  to  be  on  the  azimuth  axis.  The  sum  channel  response  indicates 
a  reduction  in  the  relative  amplitudes  of  some  of  the  targets  as  compared  to 
the  input  amplitudes.  This  effect  is  caused  by  the  beam  shape  function 
employed  in  the  simulation,  which  reduces  the  response  for  off-axis  targets. 
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Figure  18.  Effect  of  30- 
dB  system  signal-to- 
noise  ratio  on 
detei  mining  off-axis 
target  position  using 
monopulse  calibration 
on  high-range- 
resolution  target  image 
data. 


Figure  19.  Effect  of  20- 
dB  system  signal-to- 
noise  ratio  on 
determining  off-axis 
target  position  using 
monopulse  calibration 
on  high-range- 
resolution  target  image 
data. 
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The  addition  of  system  noise  has  the  same  effect  on  the  sum  channel  noise 
level  as  in  the  previous  system  noise  simulations.  For  the  off-axis  response, 
the  increase  in  noise  is  most  noticeable  for  the  targets  located  on-axis. 
Because  on-axis  targets  have  a  very  small  difference  channel  signal  level,  the 
effective  signal-to-noise  ratio  in  the  difference  channel  is  very  low  compared 
to  the  sum  channel.  This  effect  is  apparent  if  one  examines  the  baseline 
variations  of  the  off-axis  signal  for  the  on-axis  targets  as  the  system  noise  level 
is  increased.  At  a  system  signal-to-noise  ratio  of  40  dB,  the  on-axis  targets 
produce  an  almost  imperceptible  off-axis  error  signal  that  becomes  more  and 
more  detectable  as  the  signal-to-noise  ratio  is  decreased  to  20  dB. 

This  monopulse  processing  technique  can  be  quite  useful  in  resolving 
extended  targets,  because  certain  parts  of  the  target  can  be  tracked  independ¬ 
ently. 

4.3  ISAR  Results 

The  ISAR  simulations  were  done  using  the  same  target  array  as  used  in  the 
monopulse  simulations  discussed  in  section  4.2.  The  signal  returned  from  the 
rotating  target  array  was  calculated  as  discussed  in  section  3.3.  As  in  the 
monopulse  calculations,  these  simulations  were  performed  using  system 
noise  only.  Because  the  other  noise  sources  are  special  cases  of  system  noise, 
it  was  not  necessary  to  repeat  all  the  calculations  for  either  the  monopulse  or 
the  ISAR  simulations.  Separate  simulations  were  run  using  single-pulse 
signal-to-noise  ratios  of  40,  30,  20,  and  10  dB.  In  addition  a  baseline 
simulation  was  run  with  no  additive  noise. 

Figure  20  shows  the  baseline  ISAR  simulation  calculated  as  discussed  in 
section  2.2.  The  figure  shows  a  hidden-line  three-dimensional  ISAR  plot  of 
the  input  target  array.  The  range  extent  goes  from  90  to  1 05  m,  while  the  cross¬ 
range  image  covers  10.24  m.  The  amplitude  scale  goes  from  0  dBsm  to  -50 
dBsm,  with  a  total  dynamic  range  of  50  dB.  The  dynamic  range  was  limited 
by  the  inability  of  the  hidden-line  program  to  handle  multiple  overlapping 
lines.  The  target  array  was  rotated  clockwise  so  that  the  right  half  of  the  plot 
corresponds  to  positive  Doppler  frequencies  while  the  left  half  is  for  negative 
Doppler  frequencies.  The  points  that  are  off-axis  are  reduced  in  amplitude  due 
to  the  beam  shape  function  used  to  generate  the  raw  data.  With  the  dynamic 
range  limited  to  50  dB,  the  processing  side  lobes,  which  are  about  42  dB  below 
the  peak  for  a  point  target,  are  only  visible  for  the  two  points  closest  to  the 
beam  axis. 

The  positions  and  amplitudes  for  the  peak  pixels  for  each  point  are  given  in 
table  3.  The  down-range  and  cross-range  positions  of  the  targets  are  in  good 
agreement  with  the  input  positions,  and  the  target  amplitudes  indicate  the 
reduction  due  to  the  beam  shape  function. 


Range  (m) 


Figure  20.  Three-dimensional  ISAR  image  of  target  array  for  noiseless  case. 


Table  3.  Comparison  of 
Input  Target  Location 
with  the  Amplitude 
and  Position  Derived 
from  ISAR  Image 


Input  position  Output  position  Amplitude 


Range 

(m) 

Crossrangc 

(m) 

Range 

(m) 

Crossrange 

(m) 

(dBsm) 

92.93 

0.0 

93.03 

0.0 

-0.4 

101.93 

-1.0 

102.03 

-1.0 

-2.7 

94.93 

+  1.5 

95.03 

+  1.5 

-6.5 

97.43 

-2.0 

97.43 

-2.0 

-11.3 

97.43 

+2.0 

97.43 

+2.0 

-11.4 

The  effect  of  system  noise  levels  on  the  ISAR  image  quality  was  evaluated  by 
comparing  the  background  noise  generated  in  the  ISAR  image  for  the  four 
different  noise  levels.  The  background  noise  was  averaged  over  an  area  32 
pixels  on  a  side.  This  area,  covering  the  upper  right-hand  corner  of  figure  20, 
was  selected  so  that  processing  sidelobes  generated  by  each  target  were 
excluded  as  much  as  possible. 

The  results  are  shown  in  table  4.  The  table  gives  the  input  single-pulse  signal- 
to- noise  ratio,  the  average  noise  over  the  32  x  32  pixel  area,  and  the  calculated 
output  signal-to-noise  ratio  for  each  case.  The  table  indicates  an  improvement 


in  the  output  signal-to-noise  ratio  of  approximately  34  dB.  The  theoretical 
improvement  in  signal-to-noise  ratio  for  coherently  integrating  4096  separate 
samples  of  signal  plus  noise  is  36  dB.  The  results  are  in  good  agreement  with 
this  calculation;  the  slight  difference  is  probably  due  to  the  addition  of  a  small 
amount  of  sidelobe  power  from  the  individual  points. 

The  IS  AR  image  for  a  20-dB  input  signal-to-noise  ratio  is  sh-nvn  in  figure  2 1 . 
The  noise  is  apparent,  but  it  would  not  greatly  limit  the  ability  to  resolve 
targets  of  different  amplitudes  in  the  image.  Targets  as  much  as  30  dB  below 
the  peak  target  amplitude  could  still  be  distinguished  from  the  background 
noise.  Because  the  sidelobe  level  limits  the  dynamic  range  to  about  40  dB  for 
scattering  points  physically  close  together,  an  input  single-pulse  signal-to- 
noise  ratio  of  20  dB  is  probably  adequate  for  this  kind  of  image  processing. 


Table  4.  Comparison 
of  Average  Output 
Noise  Level  for  ISAR 

Input  S/N 
ratio  (dB) 

Average  noise 
(dBsm) 

Output  S/N 
ratio  (dB) 

Processing  with  Input 

oo 

-77.3 

77 

Noise  Level 

40 

-71.1 

71 

30 

-64.0 

64 

20 

-54.8 

55 

10 

—44. 1 

44 

Figure  21.  Three-dimensional  ISAR  image  of  target  array  for  20-dB  input  signal-to-noise  ratio. 
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5.  Conclusions 


We  have  developed  a  radar  simulation  program  with  the  capability  of 
generating  frequency-stepped  coherent  monopulse  signature  data  from  input 
target  arrays.  The  program  allows  one  to  calculate  the  return  from  a  three- 
dimensional  target  array  for  different  elevation  and  azimuth  angles.  The 
program  has  been  demonstrated  for  a  target  array  of  point  reflectors,  but  it 
could  be  used  to  calculate  target  signatures  from  target  model  data.  In 
addition,  system  parameters  such  as  antenna  size,  beam  squint  angle,  radar 
frequency,  frequency  step  size,  and  the  number  of  frequency  steps  per  ramp 
can  be  changed  at  will.  By  varying  these  parameters,  different  radar  systems 
can  be  simulated. 

We  have  applied  FFT  processing  techniques  to  the  simulated  frequency- 
stepped  data  to  generate  high-range-resolution  images  of  target  arrays.  These 
techniques  use  padding  to  minimize  range  response  variations  and  window¬ 
ing  to  improve  sidelobe  performance.  In  addition,  various  noise  sources  have 
been  added  to  the  simulated  data  to  determine  their  effects  on  range  images 
produced  with  FFT  processing.  For  our  application,  the  18-dB  processing 
gain  associated  with  this  technique  makes  it  possible  to  generate  reasonable 
image  quality  with  single-pulse  input  signal-to-noise  ratios  as  low  as  10  dB. 
For  applications  where  a  number  of  individual  scatterers  with  a  wide  range  of 
amplitudes  need  to  be  resolved,  a  slightly  higher  input  signal-to-noise  ratio 
would  be  required.  We  have  also  applied  these  high-range  resolution  methods 
to  the  monopulse  azimuth  channel  signal  to  allow  the  calculation  of  the  off- 
axis  position  of  scatterers  within  the  target  array.  This  technique  could  be 
applied  to  monopulse  data  collected  while  tracking  a  target,  to  select  a 
particular  part  of  the  target  to  track.  Alternatively,  this  method  could  be  used 
to  separate  the  target  sum  and  difference  signals  in  range  from  the  background 
signals.  This  would  allow  more  accurate  tracking  of  a  target  in  a  high  clutter 
environment. 

Using  the  target  simulation  program,  we  have  generated  rotating  target 
signature  data  and  calculated  two-dimensional  ISAR  images  of  the  target 
array.  We  have  developed  an  ISAR  image-processing  program  to  analyze 
these  data  and  any  other  coherent  frequency-agile  turntable  data  that  are 
available.  This  program  uses  sets  of  frequency-stepped  data  at  different 
aspect  angles  to  generate  the  images.  Since  the  input  array  consists  of  4096 
different  looks  at  the  target,  there  is  a  36-dB  improvement  in  the  signal-to- 
noise  ratio  after  processing.  This  means  that  even  at  single-pulse  signal-to- 
noise  ratios  of  20  dB,  the  image  quality  is  not  significantly  degraded.  This  was 
demonstrated  by  generating  ISAR  images  of  a  target  array  with  an  input 
signal-to-noise  ratio  of  20  dB.  The  generated  image  accurately  reproduced 
both  the  amplitude  and  position  of  the  various  scatterers  in  the  input  array. 
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We  have  developed  and  demonstrated  a  set  of  programs  for  generating  and 
processing  coherent  frequency-agile  radar  data  for  use  primarily  in  target 
signature  analysis.  These  programs  will  be  used  to  analyze  target  signature 
data  that  will  be  collected  with  the  Harry  Diamond  Laboratories  frequency- 
agile  fully  polarimetric  coherent  monopulse  95-GHz  radar. 
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Appendix  A.  Listing  of  stepped-frequency  monopulse 
signal  calculation  program 
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Appendix  A 


program  isargen64_plus 


*********************************************************** 

This  program  confutes  phases  and  anplitudes  from  * 

point  sources.  The  output  of  the  program  is  saved  * 

to  a  disk  file.  * 

Version  1.0  (NWW)  written  for  HP- 10 00  3/86  * 

Version  2.0  (siggen)  written  for  C-l  3/87  * 

Version  2.1  written  for  64x64  arrays  4/87  * 

Version  2.2  depression  and  azimuth  angles  4/87  * 
Version  2.5  added  noise, output  format  simplified 
Version  3.0  Noise  delegated  to  freq-loop  4/88  * 

Version  3.1  Inter-ramp  dead-time  added  4/88  * 

Written  by  Jerry  Silvious  ,  Branch  11300  HDL  * 

* 

Known  Bugs:  Does  not  handle  depression  <>  0  properly  * 
*********************************************************** 


implicit  none 

integer *4  j 
integer *4  stdout 


! standard  output  unit 


stdout=6 
call  isargen 
write (stdout , * ) 

write (stdout, ' ("  **  ISAR  Simulation  Terminating  **") ') 
end 

************************************************************** 

ISAR  TARGET  SIGNAL  GENERATOR  * 

This  routine  computes  phases  and  amplitudes  of  point  * 

sources  for  various  frequencies  and  various  observation  * 
angles.  The  locations  of  the  points  is  read  from  a  disk  * 
file.  The  output  is  stored  to  another  disk  file.  * 

*********** 

subroutine  isargen 


implicit  none 
complex*8  cangle 
integer *2  ntype(4) 


integer*4 
integer *4 
integer *4 
integer *4 
integer *4 
integer  *  4 
integer *4 
integer *4 
integer *4 
integer *4 
inteaer*4 


h,k, l,m 

stdin 

stdout 

stderr 

uniti 

unito 

pts 

channel 
iargc, args 
stime 
irand 


! locates  angle  of  points 

! noise  type 

! counter  variables 
! standard  in 
! standard  out 
! standard  err file 
! input  unit  for  data 
! output  unit  for  data 
! number  of  point  sources 
!beam  number  counter 
Icocmand  line  arguments 
! opsys  time 

! integer  random  number  function 
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real*4  dO (4096) ,rO (4096) ,hO (4096)  Ipoint  location  buffers 


real*4  phiO 
real*4  phi 
real*4  circlrad 


•initial  aspect  angle 
•updated  aspect  angle 
•radius  from  centroid 


real *4  depression 
real *4  rho 
real *4  azimuth 
real*4  ax 
real *4  ay 
real*4  az 

real *4  ranp_cycles 
real*4  time 
integer*4  null_ranps 
real *4  steps 
integer *4  nsteps 
real *4  deg_per_ranp 
real *4  convert 
real *4  norm 


!  depression  angle 
! distance  to  center  of  rotation 
! azimuth  angle  of  antenna 
•crossrange  coordinate  of  antenna 
!  downrange  '  '  ' 

! height  of  antenna 

•Number  of  ranps  to  be  executed 

♦time  variable 

!  ranps  for  which  there  is  no  data 
I  frequency  step  size 
Inumber  of  steps 

! rotation  rate 

!  pi/180 


real*4  nanpr 
real *4  vec_phase(3) 
real *4  vec_anp(3) 
real*4  beanupeak 
real*4  new_phase 
real*4  amp (4096) 
real*4  ajit (2) ,pjit (2) 
real*8  dbuff(6) 


! noisy  amplitude  component 
!  vector  phase 
! vector  amplitude 

I  noisy  total  phase  conponent 
Ipoint  amplitude  array 
•noise  factors 
I  real -valued  databuffer 


real*8  bfreq 
real*8  warp 
real *8  freq 
real*8  lambda 


I  beam  parameters 
I  speed  of  light  (ir/sec) 
I  updated  frequency 
I  wavelength 


real*8  pi 

real *8  theta 

real *8  hypotenuse 

real*8  new_d0 

real *8  aO ,  radius , height 

real*8  squint 

real*8  dia 

real *8  two 


Iconputed  phase  angle 
I  radial  distance  for  phase  calculation 
I  updated  downrange  location 
!  updated  point  locations 


real*4  rnd (6) .distribution (4096) 

real *4  rand 

real*4  unique_amp  (3) 

real*4  unique_phase ( 3) 

real *4  outbuf f (3072) 

character  infile*80, outfile*80 


I  random  number  collection 
lopsys  random  number  function 
I  system  noise  component 
I  random  system  noise  phase 
I  output  buffer 


parameter 

parameter 

parameter 

parameter 

parameter 


(  stdin  =  5  ) 

(  stdout=  6  ) 

(  stderr  =  0  ) 
(  uniti  =10  ) 

(  unito  =  20  ) 
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parameter 

parameter 

parameter 

parameter 

parameter 

parameter 

parameter 

parameter 


(  warp  =  299792458.  ) 

(  convert  =  .017453292) 

(  pi  =  3.14159265358979  ) 

(  bfreq  =  94.44E9  ) 

(  squint=7 .443745e-3  ) 

(  dia=. 15748  ) 

(  dia=. 0381  ) 

(  null_ramps  =  8) 


!  pi/180. 


!  transmit  frequency  95  GHz 
!  .3125*. 02382  (see  below) 
!6.2"  antenna  diameter 
!1.5"  antenna  diameter 
!  no  data  during  2  ramps 


common  /specs/hypotenuse , ax , ay, az , aO ( radius, height 


two  =  dble(2.) 


c*  *  *  *  *  Recover  command  line  arguments 


args=iargc() 
if  (args.lt. 2)  then 

write (stderr, *)  'Usage:  isargen64+  input_file  output_file' 

stop 
endif 


write (stderr, *) 

write  (stderr,  '  ("********  ISAR  SIGNAL  GENERATION  ********")') 
call  getarg(l, infile) 

open(uniti, file=infile, status='old' , READONLY) 
call  getarg(2,outfile) 

open (unito, f i le=out  file, STATUS^ ' unknown ' ) 


write (stderr, *) 

write (stderr, ' ("Radial  distance  to  target  (m) 
read(uniti, *)  rho 
write (stderr , *)  rho 

write (stderr, ' ("Depression  angle  (deg) 

read(uniti, *)  depression 

write fstderr, *)  depression 

write (stderr, ' ("Azimuth  angle  (deg) 

read(uniti, *)  azimuth 

write (stderr, *)  azimuth 

write (stderr, ' ("Rotation  speed  (deg/ramp-ramp) 
read(uniti, *)  deg_per_ranp 
write (stderr , *)  deg_per_ramp 


write (stderr ,' ("Ramp  cycles  |  ",$)') 

read(uniti,  *)  rairp_cycles 
write (stderr, *)  ranp_cycles 

wr ite (stderr, ' ("Frequency  step  size  (MHz)  |  ",$)') 

read(uniti,  *)  st^is 
write  (stderr , *)  steps 

write (stderr ,' ("Number  of  steps  per  ramp  |  ",$)') 

read (uniti , *)  nsteps 
write (stderr , *)  nsteps 
write (stderr, *) 

write (stderr, ’ ("*********************************") ') 
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write (stderr , ' ("Type  return  to  continue. . .") ') 
read (stdin, *) 

read(uniti( *)  pts  !read  the  number  of  points 

do  k=l,pts  ’read  in  the  point  coordinates 

read (uniti, *)  dO (k) , rO (k) ,h0 (k) , amp (k) 
enddo 

read  (uniti,  *)  ntype(l)  ,  ntype(2)  ,  ntype(3)  ,ntype(4) 

read (uniti, *)  ajit  (1) ,pjit (1) ,ajit(2) ,pjit (2)  ! jitter  values 

close (uniti) 

open ( 1 0, file='randccn_nos' , READONLY)  !load  random  distribution 

do  k—1,4096  ! array 

read (10,*)  distribution (k) 
enddo 
close (10) 

c  *************  begin  data  generation  **************************** 

time=-l.  *  ( (deg_per_ranp)/nsteps/(null_ramps+l) )  '.set  time  increment 

! antenna  location 

c  ax=rho*cos (depression *pi/180. )  *sin(azimuth*pi/180.)  !X-range 

c  ay=rho*cos  (depression *pi/180  . )  *cos  (azimuth*pi/180 . )  !D-range 

c  az=rho*sin  (depression*pi/180 .)  (height 

ax  =  0 
ay  =  rho 
az  =  0 

do  itfI , ramp_cycles  (number  of  angles  loop 

wr ite (stderr , ' (al6, i4, al, $) ')  ‘Computing  Ramp  #' ,m,char  (13) 

time=(m-l)  *deg_per_rarrp  !  rotation  w/  no-data 

do  l=l,nsteps  (number  of  frequency  pulses  loop 

freq=bfreq+ float ( 1 - 1)  *steps*l ,E6  !step  the  frequency 

lambda=warp/freq  ! compute  the  wavelength 

do  h=l,6 

dbuff(h)=0.  ! clear  data  buffer 

enddo 

time=time+  ( (deg_per_ramp)  /nsteps/  (null_ranps+l) ) 
c  phi=+l . *time*pi/180 .  !  CW  rotation  \aspect  angle 

phi=-l. *time*pi/180 .  !CCW  rotation  /  (?) 

do  k=l,pts  !#  of  points  loop 

c  **********  conpute  the  return  phase 
new_d0=d0 (k) 

c  new_d0=d0 (k) +(. 234212857) * (1-1) / (160 . *nsteps)  ! walking  point 

circlrad=sqrt (new_dO*new_dO-*-rO (k) *r0 (k) ) 

cangle=cnplx(dble(rO (k) ) ,new_d0)  (compute  the  initial  angle 
call  kwik (cangle,phiO) 

if  (circlrad.eq . 0 . )  phi0=0 .  (point  at  center  of  rot. 

phiO  =  phiO  ♦  azimuth*pi/180 .  !  user  rotational  offset 
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aO=circlrad*dsin (dble (phi-phiO) )  lupdate  D-range  coord. 

radius=circlrad*dcos (dble (phi-phiO) )  lupdate  v-range  coord. 

height=hO (k)  I (cylindrical  trajectory) 

hypotenuse= (ay-aO) * (ay-aO) + (ax-radius) * (ax-radius) 

&  +  (az -height) * (az -height) 

thetadsqrt (hypotenuse) * (two* (two*pi/lambda) ) 
thetadmod (theta,  (two*pi) )  !phase,mod  2*pi 

call  nominal  (freq,  airp  (k)  ,  theta,  dbuf  f)  Icoirpute  the  anplitudes 
enddo 

beanv_peak=norm  (freq,  squint,  dia)  ! compute  noise-ref  factor 

cangle=cnplx (sngl (dbuf f (1) ) , sngl (dbuf f (2) ) ) 
call  kwik (cangle, vec_phase(l) ) 

vec_anp  (l)=sngl  (dsqrt  (dbuff  (1)  *dbuf  f  (1)  +dbuff  (2)  *dbuf  f  (2) ) ) 

cangle=cnplx (sngl (dbuff (3) ) , sngl (dbuff (4) ) ) 
call  kwik (cangle , vec_phase ( 2 ) ) 

vec_airp  (2)  =sngl  (dsqrt  (dbuf  f  (3)  *dbuf  f  (3)  +dbuf  f  (4)  *dbuf  f  (4) ) ) 

cangle=cnplx (sngl (dbuff (5) )  ,  sngl (dbuff (6) ) ) 
call  kwik (cangle, vec_phase (3) ) 

vec_amp (3) =sngl (dsqrt (dbuff (5) *dbuf f (5) +dbuf f (6)  *dbuf f (6) ) ) 

c  *******************  compute  noise  factors  ****** 

rnd(l)  =  m*float (irand (stime  ())) /I .E6+ float (1-1) 
rnd  (2)  =  m* float  (irand (stime  () )) /I .E6+64 . afloat  (1-1) 
rnd(3)  =  m*float (irand (stime  ())) /I .E6+128.  +  float (1-1) 
rnd  (4)  =  m*float (irand (stime ())) /l .E6+192 . +float (1-1) 
rnd  (5)  =  m* float (irand (stime ())) /I .E6+256. + float (1-1) 
rnd  (6)  =  m* float (irand (stime  ()) )/l .E6+320. +float (1-1) 
rnd (1) distribution (int (4096 . *rand (l.E6*int (rnd (1) ) )  +  .5) ) 
rnd  (2)  distribution  (int  (4096 .  *rand  (l.E6*int  (rnd  (2) ) )  ♦  .5) ) 
rnd (3)  distribution  (int  (4096 .  *rand  (l.E6*int  (rnd  (3) ) )  +  .  5) ) 
rnd  (4)  =rand(l  .E6*int  (rnd  (4) ) ) 
rnd  (5)  =rand(l  .E6*int  (rnd  (5) ) ) 
rnd (6) =rand(l .E6*int (rnd (6) ) ) 

if  (ntype(l) .eq.l)  then  !pulse-to-pulse  anplitude 

nanpr=rnd(l)  *10*  * (-ajit (1) /20 .) 
else 

nanpr=0 . 
endif 

if  (ntype (2) .eq . 1)  then  Ipulse- to -pulse  phase 

new_phase="dble  (pjit  (1)  *md(l)  ‘convert) 
else 

new_phase=0 . 
endif 

if  (ntype  (3)  .eq.  1)  than 

unique_arrp  (1)  =rnd  (1)  *10**^-ajit(2)/20.) 
unique_aiTp  (2)  =rnd  (2)  *10**(-ajit(2)/20.) 
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unique_anp  (3)  =rnd  (3)*10**(-ajit(2)/20.) 
unique_phase  (1)  =  md(4)*2.*pi 
unique_phase  (2)  =  md(5)*2.*pi 
un ique_phase  (3)  =  md  (6)  *  2  .  *pi 
else 

unique_arrp  (1)  =0  . 
unique_anp  (2) =0 . 
unique_anp  (3) =0 . 
unique._phase  (1)  =0  . 
unique_phase  (2) =0  . 
unique^phase (3) =0 . 
endif 
******* 

dbuf  f  (1)  =  (vec_arrp  (1)  +nampr)  *cos  (vec_phase  (1)  +new_phase) 
&  +unique_amp (1) *cos (unique_phase (1) ) 

dbuf  f  (2)  =  (vec_anp  (1)  +naiTpr)  *sin  (vec_phase  (1)  +new_phase) 
&  +unique_amp  (1)  *sin(unique_phase  (1) ) 

dbuf  f  (3)  =dbuf  f  (3) 

&  +unique_amp (2) *cos (unique_phase (2) ) 

dbuf  f  (4)  =dbuf  f  (4) 

&  +unique_aiq?  (2)  *sin  (unique_phase  (2) ) 

dbuff  (5)=dbuff  (5) 

&  +unique_anp  (3) *cos (unique _phase (3) ) 

dbuff  (6)=dbuff  (6) 

&  +unique_amp (3) *sin ( un ique_phase (3) ) 


outbuff ((1-1) *6+1) =sngl (dbuff (1))  !i  sum 

outbuff ( (1-1) *6+2) =sngl (dbuff (2))  !q  sum 

outbuff  ( (1-1)  *6+3)=sngl  (dbuff  (3) )  !i  a z 

outbuff ( (1-1)  *6+4) =sngl (dbuf f  (4) )  !q  a z 

outbuff ( (1-1)  *6+5) =sngl (dbuff (5))  !i  el 

outbuff ( (1-1) *6+6) =sngl (dbuf f (6) )  !q  el 


enddo  !end  of  frequency  loop 

write (unito, *)  (outbuff (k) ,k=l,6*nsteps) 
enddo  lend  of  aspect  angle 


close (unito) 


! close  output  file 


return 

end 


subroutine  nominal (freq, mag, theta, dbuf f) 
ixrp licit  none 


integer *4  channel 
integer *4  table (4, 4) 

real*4  mag 
real*4  anpr 
real*4  chanp 
real *4  norm 
real*4  beam_peak 

real*8  dbuff(6) 
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real *8  squint 
real *8  dia 
real*8  in,q 
real *8  freq 
real*8  theta 

data  table/1,0,0,0, 

&  0,1, 0,0, 

&  0,0, 1,0, 

&  0,0, 0,1/ 

parameter  (  squint=7 .443745e-3  )  !  .3125*. 02382  (see  below) 

c  parameter  (  dia=. 15748  )  !6.2"  antenna  diameter 

parameter  (  dia=.0381  ) 

beam_peak=norm(freq, squint , dia)  Icorrpute  sum-norm  factor 

c* *********  SUMMATION  CHANNELS  ******************* 
do  channel=l, 4 

anpr=(10**  (mag/20  . ) )  *  champ  (freq,  channel) 

in=dble  (anpr)  *dcos  (theta)  /beam_peak  !  in-phase 

q  =dble(anpr)  *dsin  (theta) /beanupeak  !  q-phase 

dbuff  (l)=dbuff  O)  +in  !i  sum  channel 

dbuff(2)=dh-  ~€  k2)  +q  !q  sum  channel 

enddo 

c**********  DIFFERENCE  CHANNELS  ******************* 
do  chan.iel=l,4 

rmpr= (10  *  *  (mag/20  . ) )  *chanp  ( freq,  channel) 

in=dble  (airpr)  *dcos  (theta)  /beauupeak  » in-phase  signal 
q  =dble  (anpr)  *dsin  (theta) /beam_peak  !  q-phase  signal 

c  azimuth  i: 

dbuf f (3) =dbuf f (3) +in*table(l, channel) +in*table (3, channel) 

&  -in*table(2 , channel) -in*table (4, channel) 

c  azimuth  q: 

dbuf f (4) =dbuf f (4) +q*table fl, channel^ +q* table ( 3, channel) 

&  -q* table (2, channel) -q* table (4, channel) 

c  elevation  i: 

dbuf f (5) =dbuf f (5) + in* tab le(l, channel) +in*table(2 , channel) 

&  -in* table (3, channel) -in*table (4, channel) 

c  elevation  q: 

dbuf f (6) =dbuf f (6)  +q*table (1 , channel) +q*table (2, channel! 

&  -q*table (3, channel) -q*table (4, channel) 

enddo 

return 

end 
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Channel  amplitude  function  * 

*************************************** 

real  function  champ  (freq, channel) 
implicit  none 


real*4  ajit 
real *4  norm 
real*4  ax 
real *4  ay 
real*4  az 

real *8  hypotenuse 
real*8  aO 
real *8  radius 
real *8  height 

real *8  freq 
real* 8  squint 
real*8  dia 
real*8  c 
real*8  psi 
real *8  length 
real *8  gamma 
real*8  pi 
real *8  by,bz,bx 
real*8  anp 

integer *4  channel 

parameter  (  squint=7 .443745e-3  ) 

parameter  (  dia=. 15748  ) 
parameter  (  dia=.0381  ) 

parameter  (  pi=3 . 14159265358979  ) 
parameter  (  c=299792458.  ) 


(amplitude  jitter 

! crossrange  coordinate  of  antenna 
(downrange  '  '  ' 

(height  of  antenna 

(absolute  range  to  point 
(relative  downrange  location 
(crossrange  location 
(height  of  point  source 

(transmitted  frequency 

(half -angle  between  beams 

! antenna  diameter 

(speed  of  light 

(angle  for  amplitude  equation 

(denominator  of  dot  product 

(dot  product  angle 

(nearest  point  on  beam 
(amplitude  (not  normalized) 

(beam  channel 

!  .3125*. 02382  half -angle  between 
!  beam  centers 
(6. 2"  antenna  diamieter 
(6. 2"  antenna  diameter 


common  /specs/hypotenuse,  ax,  ay,  az ,  aO ,  radius ,  height 


by=dsqrt  (hypotenuse) 
bz=abs  (by*dtan (squint) ) 
if  (channel. eq.l)  bx=bz 
if  (channel .eq. 2 )  bx=-bz 
if  (channel .eq. 3)  then 
bx=bz 
bz=-bz 
end  if 

if  (channel .eq. 4)  then 
bz=-bz 
bx=bz 
endif 


(use  range  calculation 
!height=of f-axis  distance 
(determine  closest  point 
!  on  the  given  beam 


length-dsqrt  (bx*bx+by*by+bz*bz)  *dsqrt  (hypotenuse)  (denominator  .product 
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if  (length. eq.O)  then 
gamma- 0 . 

else 

ganma=dacos  ( (bx*  (ax-radius)  +by*  (ay-aO)  +bz*  (az-height) ) 
&  /length) !dot  product  angle 

endif 


psi=pi*dia*dsin (gamma)  /(c/freq) 

amp=dsin  (psi+pi/2  .)  /  (psi+pi/2  .)  Samplitude 
&+dsin (psi-pi/2 . ) / (psi-pi/2  . ) 

champ=(amp*pi/4.)  *  (anp*pi/4.)  Snormalized  beam  amplitude 

return 

end 


c  ***************  Normalization  Factor  ************ 
real  function  norm(freq, squint, dia) 
implicit  none 

real *8  freq  ! transmitted  frequency 

real*8  squint  Sbeam  squint  angle 

real *8  dia 

real*8  bz 

real*8  gamma, psi 

real*8  c 

real*8  pi 

real*  8  amp 

real*8  len 

parameter  (  pi  =  3.14159265358979  ) 
parameter  (  c  =  299792458  ) 

bz=10  0  0 . *dtan (dble (squint) ) 
len=1000  .  *sqrt  (bz*bz+1000  .  *1000  .  +bz*bz) 
gamma=dacos ( (1000 . *1000 .) /len) 
psi=pi*dia*sngl (dsin  (gamma) ) / (c/freq) 

amp=dsin  (psi+pi/2  . )  /  fpsi+pi/2  . )  +dsin  (psi-pi/2  . )  /  (psi-pi/2  . ) 
norm=4 .  *  (amp*pi/4 . )  *  (amp*pi/4  . ) 

return 

end 


c*  ******************************************* 
c  KWIK  * 

c  This  routine  computes  the  corrplex  * 

c  angle  of  a  given  complex  number.  * 

c* ******************************************* 
subroutine  kwik(cangle,phiO) 


complex *8  cangle 
real*4  phiO 
real*4  pi 
real*4  rO 
real*4  dO 


!  input  corrplex  number 
!  output  angle 

Sreal  part  of  cangle 
Simaginary  part  of  cangle 
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pi=3 . 14159265358979 

rO=real (cangle) 
dO=aimag (cangle) 

if  ( (rO .eq. 0 . ) .and . (dO . ge. 0 . 
phi0=pi/2 . 
goto  100 
end  if 

if  ( (rO .eq. 0 . ) .and . (dO . It . 0 . 
phi0=3*pi/2 . 
goto  100 
end  if 

phi0=atan (dO/rO) 

if  ( (rO . gt . 0 . ) .and . (dO . ge. 0 . 
if  ( (rO . It . 0 . ) .and. (dO . ge. 0 . 
if  ( (rO . It . 0  . )  .and . (dO . It . 0 . 
if  ( (rO . gt . 0 . ) .and . (dO . It . 0 . 

100  return 
end 


) )  then 


) )  then 


))  phi0=phi0 
))  phi0=phi0+pi 
) )  phi0=phi0+pi 
))  phiO=phiO+2*pi 
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COMMANDER 

US  ARMY  MATERIEL  COMMAND 
ATTN  AMCDM,  DIR  FOR  PLANS  &  ANALYSIS 
5002  EISENHOWER  AVE 
ALEXANDRIA,  VA  22333-0001 

DIRECTOR 

US  ARMY  MISSILE  LABORATORY 
ATTN  AMSMI-RPR,  REDSTONE  SCIENTIFIC 
INFO  CENTER 

ATTN  AMSMI-RPT,  TECHNICAL  INFORMATION 
DIV 

ATTN  AMSMI-RD,  SYS  SIMULATION  *  DEV 
DIR 

ATTN  AMSMI-RD,  ADVANCED  SENSORS  DIR 
ATTN  AMSMI-RD-AS-MM,  G.  EMMONS 
ATTN  AMSMI-RD-AS-RA,  J.  LOOMIS 
REDSTONE  ARSENAL,  AL  35809 

DIRECTOR 

VIGHT  VISION  &  ELECTRO-OPTICS 
CENTER 

ATTN  TECHNICAL  LIBRARY 
ATTN  CNVEO-AC,  ADVANCED  CONCEPTS 
DIV 

FT  BELVOIR,  VA  22060 
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DISTRIBUTION  (cont'd) 


US  CHIEF  ARMY  RESEARCH  OFFICE 
(DURHAM) 

ATTN  DRXRO-EG,  DIR  ENGINEERING  DIV 
PO  BOX  1 221 1 

RESEARCH  TRIANGLE  PARK,  NC  27709 


GEORGIA  TECH  RESEARCH  INSTITUTE 
GEORGIA  INSTITUTE  OF  TECHNOLOGY 
RADAR  &  INSTRUMENTATION  LABORATORY 
ATTN  T.  L.  LANE 
ATLANTA,  GA  30332 


COMMANDER 

US  ARMY  TEST  &  EVALUATION 
COMMAND 

ATTN  STEWS-TE-LG,  S.  DICKERSON 

ATTN  STEWS-TE-AF,  F.  MORENO 

WHITE  SANDS  MISSILE  RANGE,  NM  88002 

DIRECTOR 

NAVAL  RESEARCH  LABORATORY 
ATTN  2600,  TECHNICAL  INFO  DIV 
WASHINGTON,  DC  20375 

COMMANDER 

NAVAL  SURFACE  WARFARE  CENTER 
ATTN  DX-21  LIBRARY  DIV 
DAHLGREN,  VA  224148 

COMMANDER 

NAVAL  SURFACE  WARFARE  CENTER 
ATTN  E-43,  TECHNICAL  LIB 
WHITE  OAK,  MD  20910 

COMMANDER 

NAVAL  WEAPONS  CENTER 
ATTN  38,  RESEARCH  DEPT 
ATTN  381 ,  PHYSICS  DIV 
CHINA  LAKE,  CA  93555 

HEADQUARTERS 
ARMAMENT  DIV  (AFSC) 

ATTN  MSD/ENG,  MAJ  M.  TURK 
EGLIN  AFB,  FL  32542 

DIRECTOR 

NASA 

GODDARD  SPACE  FLIGHT  CENTER 
ATTN  250,  TECH  INFO  DIV 
GREENBELT,  MD  20771 

DIRECTOR 

NASA 

LANGLEY  RESEARCH  CENTER 
ATTN  TECHNICAL  LIBRARY 
HAMPTON,  VA  23665 

ENVIRONMENTAL  RESEARCH  INSTITUTE 
OF  MICHIGAN 
ATTN  I RIA  LIBRARY 
PO  BOX  61 8 

ANN  ARBOR,  MI  48107 


LINCOLN  LABORATORY 
MASSACHUSETTS  INSTITUTE  OF 
TECHNOLOGYU 
ATTN  R.  M.  BARNES 
LEXINGTON,  MA  02173 


SANDIA  NATIONAL  LABORATORIES 
PO  BOX  5800 

ALBUQUERQUE,  NM  87185 

SIMULATION  TECHNOLOGIES,  INC 
ATTN  A.  V.  SAYLOR 
PO  BOX  7009 
HUNTSVILLE,  AL  35807 

TASC 

ATTN  K.  L.  WEEKS 

907  MAR-WALT  DRIVE 

FT  WALTON  BEACH,  FL  32548 


US  ARMY  LABORATORY  COMMAND 
ATTN  TECHNICAL  DIRECTOR,  AMSLC-TD 
ATTN  DIRECTOR,  AMCLD-SE  (N.  BERG) 

INSTALLATION  SUPPORT  ACTIVITY 
ATTN  LEGAL  OFFICE,  SLCIS-CC 

USAISC 

ATTN  RECORD  COPY,  AMSLC-IM-TS 
ATTN  TECHNICAL  REPORTS  BRANCH, 
AMSLC-IM-TH  (3  COPIES) 


HARRY  DIAMOND  LABORATORIES 

ATTN  D/DIVISION  DIRECTORS 

ATTN  DIRECTOR,  SLCHD-D 

ATTN  LIBRARY,  SLCHD-TL  (3  COPIES) 

ATTN  LIBRARY,  SLCHD-TL  (WOODBRIDGE) 

ATTN  CHIEF,  SCIENTIST,  SLCHD-CS 

ATTN  DEPUTY  DIRECTOR,  SLCHD-DD 

ATTN  ASSOC  DIRECTOR,  SLCHD-PO 

ATTN  DIRECTOR,  SLCHD-ST 

ATTN  DIRECTOR,  SLCHD-NW 

ATTN  DIRECTOR,  SLCHD-TA 

ATTN  DIRECTOR,  SLCHD-TS 

ATTN  CHIEF,  SLCHD-ST-AR 

ATTN  CHIEF,  SLCHD-ST-SS 

ATTN  CHIEF,  SLCHD-ST-MW 

ATTN  CHIEF,  SLCHD-ST-OP 

ATTN  CHIEF,  SLCHD-ST-SP 
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DISTRIBUTION 

HARRY  DIAMOND  LABORATORIES 
(cont'd) 

ATTN  CHIEF,  SLCHD-ST-R 

ATTN  CHIEF,  SLCHD-ST-SA 

ATTN  CHIEF,  SLCHD-ST-MS 

ATTN  CHIEF,  SLCHD-ST-AP 

ATTN  R,  BERG,  SLCHD-NW-EH 

ATTN  S.  KHAN,  SLCHD-NW-EH 

ATTN  W.  J.  PATTERSON,  SLCHD-NW-EH 

ATTN  P.  WALKER,  SLCHD-ST 

ATTN  D.  R.  COOK,  SLCHD-ST-AR 

ATTN  D.  L.  RODKEY,  SLCHD-ST-MS 

ATTN  J.  DAVID,  SLCHD-ST-R 

ATTN  C.  ARSEM,  SLCHD-ST-SA 

ATTN  J.  COSTANZA,  SLCHD-ST-SA 


(cont'd) 

HARRY  DIAMOND  LABORATORIES 
(cont'd) 

ATTN  D.  M.  HULL,  SLCHD-ST-SA 

ATTN  H.  LE,  SLCHD-ST-SA 

ATTN  J.  NEMARICH,  SLCHD-ST-SP 

ATTN  J.  DROPKIN,  SLCHD-ST-SP 

ATTN  G.  GOLDMAN,  SLCHD-ST-SP 

ATTN  D.  WIKNER,  SLCHD-ST-SP 

ATTN  D.  HUTCHINGS,  SLCHD-ST-SP 

ATTN  M.  PATTERSON,  SLCHD-ST-SS 

ATTN  J.  SILVERSTEIN,  SLCHD-ST-SP 

ATTN  R.  KAPOOR,  SLCHD-TA-AS  (2  COPIES) 

ATTN  R.  D.  GOODMAN,  SLCHD-TA-ES 

ATTN  J.  SILVIOUS,  SLCHD-ST-SP  (20  COPIES) 

ATTN  R.  WELLMAN,  SLCHD-ST-SP  (20  COPIES) 
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